jAs// 7 >/- //sr? 


NASA-TM-87599 19860007903 

NASA Technical Memorandum 87599 


A CONTROL SYSTEM DESIGN APPROACH 
FOR FLEXIBLE SPACECRAFT 


Larry M. Silverberg 


November 1985 


FOB 




/-at? 




^SaiSS B.OOS* 




7 


JAN 2 7 


LANGLcY RESEARGi I -KM • f .U 
LIBRARY, NASA 
HAMPTON, VIRGINIA 

NASA 

National Aeronautics and 
Space Administration 

Langley Research Center 

Hampton, Virginia 23665 



Abstract 


A control system design approach for flexible spacecraft is presented. 

The control system design is carried out in two steps. The first step 
consists of determining the "ideal" control system in terms of a desirable 
dynamic performance. The second step consists of designing a control system 
using a limited number of actuators that possess a dynamic performance that is 
close to the ideal dynamic performance. The effects of using a limited number 
of actuators is that the actual closed-loop eigenvalues differ from the ideal 
closed-loop eigenvalues. A method is presented to approximate the actual 
closed-loop eigenvalues so that the calculation of the actual closed-loop 
eigenvalues can be avoided. Depending on the application, it also may be 
desirable to apply the control forces as impulses. The effect of digitizing 
the control to produce the appropriate impulses is also examined. 
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Introduction 

A great deal of work has been dedicated to the development of structural 
control theories. Although the developments are extensive, the designer often 
finds it difficult to apply many of these theories to "real structural control 
problems." Indeed, it is of primary concern to bridge the gap between 
engineering design and the structural control theories. 

The focus of the research into structural control theories is diverse. 
Many of the researchers are concerned with practical implementation problems 
and toward that end, promote decentralized control (Refs. 1 and 2). Others 
promote a centralized modal control approach and toward that end, point out 
that a control theory should not destroy certain characteristics which are 
natural to a structure (Refs. 3-5). Still others, in search for a global 
optimum, are concerned with distributed controls (Refs. 6 and 7). Much 
attention is also given to describing the robustness of the control theories 
in the presence of modelling errors, particularly in view of the fact that it 
is difficult to characterize structural stiffness in mathematical models 
(Refs. 8-10). 

All of these concerns support the objective to uniformly dampen the 
motion of a spacecraft. As it turns out, a uniform damping control is a 
robust, decentralized, natural control with near globally optimal performance 
(Ref. 11) . Thus, a uniform damping control answers the concerns raised in the 
previously cited references. 

In this paper, an engineering design approach to structural control is 
described. The design of a uniform damping control system is carried out in 
two independent steps. The first step consists of identifying the solution 
which leads to the ideal dynamic performance. Toward that end, one recognizes 
that the state of a spacecraft is distributed over its domain, implying that 
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the ideal dynamic performance will require distributed actuation and sensing 
devices. On the other hand, it is recognized that the use of these 
distributed devices is, for the most part, impractical. The second step 
consists of constructing a control system of minimal cost which exhibits 
dynamic performance that is as close as possible to the ideal. Therefore, the 
second step consists of implementing the uniform damping control obtained in 
the first step using discrete actuation and discrete sensing devices. As it 
turns out, ideal performances can be obtained with a relatively small number 
of actuators. 

II. Mathematical Description 

The equations of motion of a flexible structure can be expressed in the 

form 

Mx(t) + Kx(t) = F(t) (1) 

where x(t) is an n-dimensional vector of nodal displacements and slopes and 
F(t) are forces and moments at the corresponding nodes. M and K denote n by 
n mass and stiffness matrices, respectively, and overdots represent 
differentiations with respect to time. The mass and stiffness matrices are 
obtained using the finite element method. Common computer programs capable of 
generating the mass and stiffness matrices include NASTRAN and SAP. 

Associated with the equations of motion, one commonly defines the 
eigenvalue problem 

= Kj, (2) 

The solution of this problem is known as the eigensolution which consists of 
the eigenvector £ and the associated eigenvalue x. There exist n 
eigensolutions, i.e. n eigenvectors £ p (r = 1, 2, ..., n) and n associated 
eigenvalues x p (r = 1, 2, ..., n). Structural dynamicists commonly refer to 
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the eigenvectors as natural modes of vibration . The associated eigenvalues 
are related to the natural frequencies w p by x p = w p 2 (r = 1, 2, ..., n). As 
a general rule of thumb, the computed eigensolution with higher associated 
natural frequencies are inexact. Indeed, modelling error will significantly 
effect these quantities. Only the eigensolutions with lower associated 
natural frequencies can be computed accurately. However, more often than not, 
we are only concerned with the lower modes, so this presents no difficulty. 
NASTRAN and SAP are two typical computer programs capable of computing the 
eigensolution (Ref. 12). 

The natural modes can be normalized so that 

*X* Vs < 3 > 

where S ps = 0 for r 4 s and 6 pp = 1. We express the displacement vector X(t) 
as a linear combination of the lowest m modes, written 

x(t) = ijUjft) + $ 2 (t)u 2 (t) + ... + <yj m (t) (4) 
where m<<n, and u p (t) (r = 1, 2, ..., m) are modal displacements which express 
the degree to which the modes participate in the system response. Generally, 
the higher modes do not contribute significantly in the response so they are 
not included in Equation (4). The modal displacements are governed by the 
scalar equations, 

u p (t) +u)2 Up (t) = f p (t), (r = 1, 2, ..., m) (5) 

where the modal forces f p (t) are related to the nodal force F(t) by 

f p (t) = F(t), (r = 1, 2, ..., m) (6) 

We have assumed here that the modes are normalized, i.e. that Equation (3) is 
satisfied. It remains to compute the modal displacements in Equation (6). 
Toward this end, we first distinguish between rigid-body modes for which 
u p = 0 and flexible-body modes for which 0. 
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(A) Rigid-body Modal Responses (w p = 0) 

We rewrite Equation (5) in the state space by introducing the change of 
variables u p (t) = [u p (t) u (t)]^ and obtain the modal equations 

u r (t) = Au p ( t ) + Bf p (t) (7) 

where 

a = [; s] • s - u ( 8 > 

The solution to Equation (8) can be converted into a difference equation. 
Letting T denote the time step, and letting u (k) and u p (k) denote the modal 
velocity and modal dispacement at time kT, (k =0, 1, 2, . ..) we obtain the 
difference equations 

u p (k + 1) = u p (k) + Tf p (k) (9a) 

u p (k + 1) = u p (k)T + u p (k) + ^f p (k) (9b) 

Equation (9) is used to compute the response of a rigid-body mode. 

(B) Flexible-body Modal Responses (w i- 0) 

Equation (5) describes the motion of an undamped oscillator. However, 
structures experience small degrees of structural damping. We can introduce 
some damping into the mathematical model at the exponential rate a p by 
replacing Equation (5) with 

u p (t) + 2a p u p (t) + («2 + uj2)u r (t) = f p (t) (10) 

The natural frequency in Equation (10) is identical to that in Equation (5). 

We rewrite Equation (10) by introducing the change of complex variables 
u p (t) = Re{w p (t)}, u p (t) = Re{A r w p (t)} where A p = -a + iu , and we obtain the 
complex modal state equations 

w p (t) = X p w p (t) + f p (t)/(iu p ) 


(ID 
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Letting T denote the time step, the response to Equation (11) is given by the 
difference equation 

w p (k+l) ■= $ r w r (k) + r p f p (k) (r = 1, 2, m) (12) 

where 

V 

* r = e ,r p = (* p - l)/(ix r u p ) (13) 

Equation (12) is used in order to compute the response of a flexible-body 
mode. For these purposes, it is desirable to take a time step smaller than 
one tenth of the smallest flexible body period of oscillation. 

III. Control System Design 

The control system design is carried out in two steps. In the first 
step, one constructs the "ideal" control system with the best dynamic 
performance that nature will allow. Such a system requires distributed forces 
which are certainly impractical for most applications. The second step 
consists of designing a control system of minimal cost and greatest simplicity 
and one which imitates the ideal control system. Perhaps the simplest way to 
carry out the second step is to consider various designs and to compare the 
dynamic performances of these designs with the dynamic performance of the 
ideal control system. 

Step 1: The ideal control system. 

For vibration suppression, pointing, and shape control, the ideal control 
system is one which dampens all the modes of vibration at a single exponential 
rate a (Ref. 11). The linear feedback control law is 

F(t) = -2aMx(t) - a 2 Mx(t) (14) 

Substituting Equation (4) into Equation (14) while considering the 
orthonormality conditions, Equation (3), we obtain the expressions for the 
modal control forces 

f (t) = -2ail r (t) - a 2 U r (t), (r = 1, 2, 


m) 


(15) 
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We observe from Equation (15) that only the rth modal displacement and the rth 
modal velocity control the rth modal force. Such a control is referred to as 
natural because the modal coordinates do not couple the equations of motion 
(Refs. 7 and 11). Substituting Equation (15) into Equation (5), we obtain the 
closed-loop modal equations 

Up ( t ) + 2aUp(t) + (a 2 + a) r 2 ) u p ( t ) = 0 (r = 1, 2 m) (16) 

The corresponding closed-loop eigenvalues are given by 

Aj p = l/2[-2a + / (2a) 2 - 4 (a 2 + ay 2 )] = -a + i ay (17) 

From Equation (17), the closed-loop modes all decay at the same exponential 
rate a and the closed-loop frequencies of oscillation are identical to the 

i 

natural frequencies. Also, observe that the control law. Equation (14), is 
independent of the spacecraft stiffness. As a general rule of thumb, when a 
control system is designed to dampen modes in a more non-uniform manner, the 
control law will tend to depend more on the structural stiffness . Therefore, 
in the interest of designing a robust control system and one which does not 
depend explicitly on the fidelity of the mathematical model of stiffness, we 
uniformly dampen the motion. 

The objective to uniformly dampen the motion can also be arrived at from 
other points of view. For example, let us assume that we wish to drive the 
motion of a given point on the structure to equilibrium at the exponential 
rate a, i.e. we wish that a given point be exponentially stable. Then, it can 
be shown that this point will be exponentially stable at the exponential decay 
rate a only if all of the natural modes of vibration are exponentially stable 
at the rates a p not less than a. Also, note that any ef-fort to dampen a given 
mode at an exponential rate a p strictly greater than a will require unnecessary 
fuel . Therefore, the most effective way to drive the motion of any point to 
equilibrium at the exponential decay rate g is by damping the motion of the 


natural modes uniformly at the exponential decay rate a (Ref. 11). 
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Finally, we observe that the uniform damping control law, Equation (14), 
is decentralized. Because the mass matrix is diagonal, if we write 

M = diag(m 1 , m 2 , m n ), (18) 

then Equation (14) becomes 

F p (t) = -2am r x r (t) - a 2 m r x p (t), (r = 1, 2 n) (19) 

Clearly, Equation (19) represents a set of independent control laws, which 
suggests that uniform damping is relatively easy to implement. As a matter of 
theoretical interest, uniform damping control represents a close approximation 
to globally optimal control (Ref. 11). 

In view of the considerations presented in the previous paragraphs, the 
objective to uniformly dampen the motion has been chosen, and for the purpose 
of design, it will be viewed herein as an ideal. 


Step 2: Implementation of the ideal control system. 

It is usually impractical to consider a large number of control forces as 
in Equation (19). Therefore, we arrive at the second step and design a 
control system that performs as closely as possible to the ideal control 
system. The control law obtained in the second step can be given by 

F(t) = -Cx(t) - Dx(t) (20) 

where C and D are usually sparce matrices because in most applications only a 
relatively small number of control forces are required. It is of immediate 
concern to describe the degradation in performance due to implementing the 
controls with a limited number of control forces. As it turns out, the 
degradation in performance can be marginal. Substituting Equation (20) into 
Equation (6) and considering Equations (3) and (4), we obtain the modal 
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Equation (21) can be rewritten in the form 

" m T . 

u p (t) + 2au p (t) + (a 2 + a> 2 )u r (t) = 2a6 ps )u s (t) + 

- a2<S rs ^ U S ( r = lj 2 ’ ( 22 ) 

The flexible-body modes and the rigid-body modes in Equation (22) can we 

rewritten in the state space by introducing the complex change of variables 

u p (t) = Re (w p (t)}, ii p (t) = Re {A p w r (t)}, (r = 1, 2, ..., m) (23) 

where A p = -a + i a> r are the system eigenvalues that would be obtained using 

the ideal control system. We obtain the complex modal state equations 

m _____ 

w r (t) = X p w (t) + 1/ 2 e (g ps w s (t) + 9 rs w (t)) (24) 

s=l 

where 


’rs 


= ( a26 rs“ ) / ( 1 ) + ( 2a6 r s" 4r C is^\^ iu r^ 


(25) 


(r,s 1, 2, • • • j m) 

The eigenvalues of the controlled spacecraft lie in the circles with centers 
C p and associated radii R p , given by 


m 

X + g rr /2, R r = S 
' ■ S=1 


TS 


(26) 


s*r 

Note that the centers C p are also first-order approximations of the 
eigenvalues associated with the ideal control system. Equation (26) can be 
used in order to compare the performance of the control system design with the 
performance of the ideal control system. 


IV. Digitization of the Controls 

In the previous section, distributed controls were discretized in space 
leading to the implementation of the controls using a limited number of 
control forces. The controls acted continuously in time. The controls can 
also be discretized in time leading to digital controls. In the process, the 
dynamic performance of the controls are expected to change depending on the 
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level of digitization. The question arises, at what level of digitization 
will the dynamic performance of the spacecraft vary significantly from the 
dynamic performance of the spacecraft with an ideal control system. Consider 
the continuous controls acting at the rth node with the associated control law 
F p (t) = -2am r x r (t)-a 2 m r x r (t) (27) 

Here, m p refers to the mass of the region within which the control force F p (t) 

acts. Over a small time increment T, we apply an impulse 

I p (t) = /t +TF r ( T ) dT = F r «0T (28) 

so that 

I p (t) = -2aTm r x r (t)-a 2 Tm r x r (t) (29) 

Instead of applying continuously acting controls as suggested by 

Equation (27), let us apply an impulse every k seconds. Then, we replace the 

continuous control law. Equation (27), with the digital control law 

I p (t) = -2aKTm r x r (t)-a 2 KTm r x r (t) (30) 

where the impulse I r (t) is applied every KT seconds. The particular effects 
of implementing Equation (30) rather than Equation (27) are described in the 
numerical example. 

V. Uniform Damping of a Simply Supported Beam 

As an illustrative example, we control a simply supported beam of length 
a = 10.0 units with unit mass per unit length and unit stiffness density. For 
this simple example, the equations of motion admit closed-form expressions* 

The normalized eigenfunctions and natural frequencies are 

* r (x' - (2/a) 1 / 2 Sln(£&) „ r = (^) 2 , r - 1, 2 m (31) 

For the sake of this example, we assume that the lowest m = 10 modes of 
vibration contribute significantly to the overall system response and that the 
contribution of the remaining modes to the motion is negligible. The beam is 
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given an initial unit step input at x = 4.0 for 2.0 seconds. We design for a 
uniform exponential decay rate of a = 1.0 and we assume that 1 percent 
structural damping is present in the beam. 

As a first step, the ideal control system is designed. The free response 
is shown in Figure 1 and the ideal control system response is shown in 
Figure 2. The ideal closed-loop eigenvalues are given in Table 1. Next we 
consider implementing the control system using a discrete number of control 
forces. In order to approximate the ideal control system, we locate control 
forces along the beam at the points P , (r = 1, 2, ..., s; s = 4, 5) 

(See Table 2). The associated control laws are given by 

F p (t) = -2am r x r (t)-oc 2 m r x r (t), m r = a/s, (r = 1, 2, ..., s) (32) 
where x p (t) is the displacement at P p . Here, again, m p represents the mass in 
the region of the rth control force. The responses of the beam with the 
discrete controls are shown in Figures 3 and 4. The corresponding fuels 
consumed by the controls are shown in Figures 5 and 6. Also, the 
corresponding first-order approximations of the closed-loop eigenvalues are 
given in Tables 3 and 4. 

Next we digitize the control law Equation (30). The responses of the beam 
using digitized discrete controls are shown in Figures 7 and 8. The 
corresponding fuels consumed by the controls are shown in Figures 9 and 10. A 
computer program listing is given in Appendix A. 

VI. Conclusions 

A control system design approach for flexible spacecraft has been 
presented. The control system design is carried out in two steps. The first 
step consists of determining an "ideal" uniform exponential rate at which we 
desire the spacecraft motion to dampen. Next, we construct a control with 



dynamic performance that is close to the "ideal" using a limited number of 
actuators. It is also shown that the controls can be digitized when it is 
desirable to create forces using impulses. 

The control system design approach is demonstrated with a simple 
numerical example in which it is shown that close to ideal dynamic 
performances can be obtained with a relatively small number of actuators. 
Also, the effects of digitizing the controls on the dynamic performance is 


illustrated. 
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iop Eigenvalues 

a + iw 

r 

-a+ iw 

r 

-1.0 + i0.098 
-1.0 + iO. 394 
-1.0 + i0.888 
-1.0 + il .579 
-1.0 + i2 . 467 
-1.0 + i3. 553 
-1.0 + i4.836 
-1.0 + 16.316 
-1.0 + 17.994 
-1.0 + 19.869 


1 



Locations P r of the Control Forces 



P 1 

P 2 

P 3 

P 4 

P 5 

Five Forces 

1.0 

3.0 

5.0 

7.0 

9.0 

Four Forces 

2.0 

4.0 

6.0 

8.0 

— 


Table 2 



First-Order Approximation of the 
Closed-Loop Eigenvalues Using Five Control Forces 


r 

X = 
r 

-a + ia> 
r r 

1 

-1.0 

+ 10.098 

2 

-1.0 

+ 10.394 

3 

-1.0 

+ 10.888 

4 

-1.0 

+ il.579 

5 

-2.0 

+ i2 .264 

6 

-1.0 

+ 13.553 

7 

-1.0 

+ i4.836 

8 

-1.0 

+ i6.316 

9 

-1.0 

+ i7.994 

10 

0.0 

+ 19.920 


Table 3 



First Order Approximation of the 
Closed Loop Eigenvalues Using Four Control Forces 


r 

X = 
r 

-a r + iw r 

1 

-1.25 

+ 10.078 

2 

-1.25 

+ i0.747 

3 

-1.25 

+ 11.167 

4 

-1.25 

+ il.500 

3 

-1.25 

+ i2.760 

6 

0.00 

+ 13.517 

7 

-1.25 

+ 14.810 

8 

-1.25 

+ i6.296 

9 

-1.25 

+ 17.978 

10 

0.00 

+ 19.920 


Table 4 



(3‘6) n 























Figure 5. Fuel-Five Control 
Forces, Continuous in Time. 


Figure 6. Fuel-Four Control 
Forces, Continuous in Time. 
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Figure 8. Controlled Response-Four Control Forces, Impulses Every 

0.3 Seconds. 









Figure 9. Fuel-Four Control 
Forces, Impulses Every 
0.2 Seconds. 


Figure 10. Fuel-Four Control 
Forces, Impulses Every 0.3 
Seconds. 





Appendix A. Computer Program Listing. 


♦NULL* 

occc ccccccccccccccccccccooc 
cccc ccccocccccccccccccccccc 
occc 

OCCC DATA GENERATION PROGRAM 
CCCC 

CCCC THE STRUCTURAL SYSTEM PARAMETERS 

CCCC ARE DEFINED INCLUDING THE NATURAL 

CCCC FREQUENCIES AND THE NATURAL MODES* 

CCCC 

CCCC CCCCCCCCCCCCCCCCCCCCCCC 
CCCC CCCCCCCCCCCCCCCCCCCCCCC 

cccc 

REALX8 UEC(25) 

t amha 

OPENCUNIT=ll,FILE='DAT',STATUS='UNKNOUN' ) 
M = 10 
N = 9 

TSTEP=0.05 
NSTEP=200 
URITEC11,X)M 
URITEC11,*)N 
UR I TEC 11,X )TSTEP 
UR I TEC 11 # X )NSTEP 
PI=ACOS(-l ♦ ) 

AA=10* 

ZETA=0*01 

SQ20A = SQRTC2./AA ) 

DO 1 1 = 1, M 

OMEGA= C IXPI/AA )XX2 



Appendix A. Continued. 


OMEGA* (I*PI/AA)**2 
ALFA=£**ZETA*OMEGA 
LAMDA= (0.^1. )*OMEGA-ALFA 
URITEC 11 ,* )LAMDA 
i CONTINUE 

DO 3 1=1, N 
DO £ J=1,M 

N1 =N+1 

UEC(J )=SG£OA*SIN(J#PI*( 1-0*0 )/Nl ) 
£ CONTINUE 

URITE( 11, 180 ) ( UEC( J ), J = l, M ) 

3 CONTINUE 

100 F0RMAT(£X,5E15*6) 

CLOSEC 11 ) 

STOP 

END 


BOTTOM 



Appendix A. Continued. 


♦NULL. 

CCCC 

cccc cccccccccccccccccccccc 
CCCC cccccccccccccccccccccc 
cccc 

CCCC EXTERNAL DISTURBANCE PROGRAM 
CCCC 

CCCC THE EXTERNAL FORCES NOT INCLUDING CONTROL 

CCCC FORCES ARE DEFINED. 

CCCC 

CCCC CCCCCCCCCCCCCCCCCCCCCC 
CCCC CCCCCCCCCCCCCCCCCCCCCC 
CCCC 

OPEN (UN I T= 13, FILE* 'FORCES' , STATUS* 'UNKNOWN' ) 

NP=1 

IFOR1 =4 

WRITE(13,X)NP 

WRITE(13,X)IF0R1 

Fl = l 

F4 = 0 

DO 1 K=l,40 

1 URITE ( 13,X )F1 
DO E K=21,200 

2 WRITE(13,X)F4 
CLOSEC 13 ) 

STOP 

END 


BOTTOM 



Appendix A. Continued. 


♦NULL* 

COCO 

COCO cccccoccoccccccccccccccc 
cccc cccccocccccccoccoccccccc 
cocc 

COCO THE CONTROL PARAMETERS ARE DEFINED 
CCCC 

CCCC cccccccccccccccccccccccc 
cccc cccccccccccccccccccccccc 
cccc 

REAL*8 XMASSC9) 

OPEN(UNIT» 11, FILE* 'CONTROL', STATUS- 'UNKNOWN' ) 
OPEN (UNXT=1 2, FILE= 'DAT' , STATUS* 'UNKNOWN' ) 
ALFA=1 *0 
READ (12,*)™ 

READ(12,*)N 
DO 18 J = i,N 
XMASSC J ) = 10*/N 
10 CONTINUE 
KTIME= 1 

UR I TEC 11,* ) (XMASS( I ), I = 1,N ) 

URITEC 11,* )KTIME 
URITEC 11 ,* )ALFA 
CLOSE ( 11 ) 

CLOSEC 12 ) 

STOP 

END 


BOTTOM 


Appendix A. Continued. 


♦NULL* 

SUBROUTINE LAU ( FOR, X, XDOT, I , FORT, XNASS, KTINE, ALFA ) 

CCCC 

CCCC CCCCCCCCCCCCCCCCCCCCCCCC 
CCCC CCCCCCCCCCCCCCCCCCCCCCCC 
CCCC 

CCCC SUBROUTINE LAU 
CCCC 

CCCC THE CONTROL LAU IS DEFINED* 

CCCC 

CCCC CCCCCCCCCCCCCCCCCCCCCCCC 
CCCC CCCCCCCCCCCCCCCCCCCCCCCC 
CCCC 

REALX8 X ( 9 ) , XDOT ( 9 ) , FOR ( 9 ) , XPIASS ( 9 ) , FORT ( 9 ) 

DO 1 K=l,9 

I F ( ( I/KTIP1E )XKTIf1E.NE. I )GOTOi 

F0RK = -ALFAXXf«lASS(K)*KTinEX(2**XD0TCK)+ALFA*X(K)) 

FOR(K )=FOR(K )+FORK 

FORT ( K )= FORT ( K ) +ABS ( FORK ) 

1 CONTINUE 
RETURN 
END 


BOTTOP1 



Appendix A. Continued. 


♦NULL. 

SUBROUTINE RESP ( UEC ,UAL,X, XDOT, T, FOR, N, N, U, UDOT, U ) 

CCCC 

ccoc cccccoccccccccoccccc 
CCCC cccccccccccccccccccc 

CCCC 

CCCC SUBROUTINE RESP 
CCCC 

CCCC THE SYSTEN RESPONSE IS UPDATED FOR EACH 

CCCC TINE STEP. THE CONPUTATION DISTIGUISHES 

CCCC BETUEEN RIGID-BODY NOTION AND FLEXIBLE-BODY 

CCCC NOTION. 

CCCC 

CCCC CCCCCCCCCCCCCCCCCCCC 
CCCC CCCCCCCCCCCCCCCCCCCC 
CCCC 

REALX8 UEC(9,25),X(9),XD0T(9 ),F0R(9),U(25 ), UDOT (25 ) 
C0NPLEXX16 UAL (25 ),U(25 ),PSI, GANA, ONI 
DO 3 J = 1,N 
F = 0 

DO 1 K = 1,N 

i F=F+UEC(K, J )XFOR(K ) 

IF(CDABS(UAL(J) ).LT.1.D-6)G0T0 2 
PSI=CDEXP(UAL(J)*T) 

ON= (0. ,-l . )XUAL(J ) 

ONI= (0. , i ♦ )XON 
GANA= (PSI-1 )/UAL( J)/ONI 
U(J)=PSIXU(J)+GANA*F 
U( J)=U(J) 

UDOT ( J ) = UAL ( J )XU ( J ) 
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UDOT(J)=UAL(J)*U( J) 

GOTO 3 

2 U( J )=U(J)+T*UD0T( J )+T#*2/a«#F 
UD0T( J ) =UD0T ( J )+TXF 

3 CONTINUE 
DO 4 K=1,N 
X ( K ) =0 
XDOT ( K ) =0 
DO 4 J*i,N 

X(K )=X(K )+UEC(K> J )*U( J ) 
XD0T(K)=XD0T(K )+UEC(K # J )*UD0T( J) 

4 CONTINUE 
RETURN 
END 


BOTTOM 
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.NULL* 

COCO 

cccc cccccocccccccccccc 
cccc occccccccccccooccc 
cccc 

CCCC SYSTEM RESPONSE PROGRAM 
CCCC 

CCCC THE RESPONSE OF THE CONTROLLED SYSTEM IS 
CCCC COMPUTED AT UARIOUS POINTS. 

CCCC 

CCCC CCCCCOCCCCCCCCCCCC 
CCCC CCCCCOCCCCCCCCCCCC 

cccc 

REALX8 UEC(9,25),X(9),XD0T(9),F0R(9),U(25),UD0T(25),F0RT(9) 
REALX8 XMASS ( 9 ) 

C0MPLEXX16 UAL(25),U(2S) 

INTEGER IF0R(9) 

OPEN(UNIT=ll ,FILE= 'DAT' , STATUS® 'UNKNOUN' ) 

OPEN (UNIT =13, FILE* 'FORCES ' , STATUS* 'UNKNOWN' ) 

0PEN(UNIT=14, FILE* 'OUTl', STATUS* 'UNKNOUN' ) 

OPEN (UNIT = 15, FILE* '0UT2' , STATUS* 'UNKNOUN' ) 

OPEN (UNIT* 16, FILE* '0UT3' , STATUS* 'UNKNOUN' ) 

OPEN (UNIT =17,FILE='0UT4' , STATUS* 'UNKNOUN' ) 

OPEN (UNIT* 18, FILE* ' OUTS' , STATUS* 'UNKNOUN' ) 

0PEN(UNIT=6, FILE* 'FORI', STATUS*' UNKNOWN 7 ) 

OPEN (UNIT =7, FI LE= 'F0R2' , STATUS = 'UNKNOUN' ) 

OPEN (UNIT = 8,FILE= 'FOR3' , STATUS* 'UNKNOWN' ) 

OPEN (UNIT = 9, FILE* ' F0R4 ' , STATUS* 'UNKNOWN' ) 

OPEN (UNIT =10, FILE* 'F0R5' , STATUS* ' UNKNOUN ' ) 

0PEN(UNIT = 19, FILE* 'CONTROL', STATUS* 'UNKNOUN' ) 
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OPEN (UNIT= 19, FILE 3 'CONTROL' , STATUS 3 'UNKNOWN' ) 

READ(11,X)N 

READ(11,X)N 

READ(11,X)T 

READ( 11,X )L 

READ(il,X)(UAL(I),I 3 l,M) 

READ(11,X)( (UEC(I,J), J = i,M) , I 3 1,N ) 
READ(19.X)(XNASS(I),I 3 1,N) 

READ<19,X)KTIf1E 
READ( 19, X )ALFA 
URITE( 14, X )L,l,L 
URITE(15,X )L,L,L 
URITE( 16, X )L,L,L 
URITE( 17, X )L,L,L 
URITEC 18,X )L,L,L 
URITE(6, X )L,L, L 
URITE(7,X)L,L,L 
URITE(8,X )L,L,L 
URITE(9,X )L,L,L 
URITE(10,X)L,L,L 
TN = 0 

DO 1 1 = 1, M 
U ( I ) = 0 
UDOT ( I ) = 0 
U(I ) = 0 
CONTINUE 
DO 2 K= 1 , N 
FORT ( K ) =0 
IFORCK ) = 0 
X(K )=0 
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# 


to 


X(K)=0 
XD0T(K ) = 0 

2 CONTINUE 
READ( 13, % )NP 

READ ( 1 3 . * ) ( I FOR ( K ) , K = 1 , NP ) 

DO 4 I = 1,L 
DO 3 fC-1 , N 
FOR(K)=0 

3 CONTINUE 

READ ( 13 , % ) (F0R( IFOR(K) ),K=i,NP ) 

CALL LAIJ ( FOR, X, XDOT, I , FORT, XMASS, KTIME, ALFA ) 
CALL RESP(UEC,UAL,X,XDOT,T,FOR,M,N,U,UDOT,U ) 
URITE( 14, 100 )TM,X( 1 ),XDOT( 1 ) 

URITE( 15, 100 )Tf],X<3 ),XD0T(3 ) 
URITE(16,100)TM,X(5),XDOT(5) 
URITE(17,100)Tn,X(7 ), XD0T(7 ) 
URITE(1S,100)TH,X(9 ),XD0T(9 ) 

URITEC6, 100 )TM,FORT ( 1 ),F0RT(2 ) 

URITE( 7, 100 )TN, FORT (3 ),FORT ( 4 ) 

URITEC8, 100 )TM,F0RT(5 ),F0RT(6 ) 

URITE(9, 100 )TM,FORT (7 ),F0RT(8 ) 

URITEC 10, 100 )TM,FGRT(9 ) ,F0RT(9 ) 

TM=T+TM 

4 CONTINUE 

100 FORMAT (F6,3,2E22* 13 ) 

CLOSEdl ) 

CL0SEC13) 

CL0SEC14) 

CL0SEC15) 

CL0SEC16) 
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CL0SE( 16 ) 
CLOSE (17) 
CLOSE (18) 
CLOSE (6 ) 
CLOSE (7) 
CLOSE (8) 
CLOSE ( 9 ) 
CLOSE (10) 
STOP 
END 


BOTTOM 
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♦NULL. 

CCCC CCCCCCCCCCCCCCCC 

occc cccccccccoccoccc 

CCCC 

CCCC CONTROL ROBUSTNESS PROGRAM 
CCCC 

CCCC IDEALLY, A DESIRABLE DYNAMIC PERFORMANCE 

CCCC REQUIRES DISTRIBUTED SENSING AND ACTUATION 

CCCC UHICH IS FOR THE MOST PART IMPRACTICAL. 

CCCC THEREFORE, ONE RESORTS TO FINITE-DIMENSIONAL 
CCCC SENSING AND ACTUATION. THIS PROCESS OF GOING FROM 
CCCC DISTRIBUTED TO DISCRETE IS CALLED CONTROL 
CCCC DISCRETIZATION. THIS PROGRAM LOOKS AT THE EFFECTS 
CCCC OF CONTROL DISCRETIZATION ON THE DYNAMIC PERFORMANCE. 
CCCC TOUARD THIS END, UE LOOK ATi 

CCCC 1) CHANGES IN THE NEIGHBOURHOODS OF THE CLOSED-LOOP 
CCCC EIGENUALUES. 

CCCC 2) FIRST-ORDER PERTUBATIONS OF THE CLOSED-LOOP 
CCCC EIGENUALUES. 

CCCC 

CCCC 

CCCC CCCCCCCCCCCCCCCC 
CCCC CCCCCCCCCCCCCCCC 
CCCC 

REALX2 UEC(25,25 ),C(25,25 ),D(25,25 ),XMASS(25 ),RAD(25 ) 
C0MPLEXX16 UAL (25 ),CEN(25 ),LAM(25 ),GRS,GJS,GJI,0M(25 ) 
OPEN ( UNIT ■ = 1 1 , F ILE* ' DAT ' , STATUS = ' UNKNOWN ' ) 

OPEN(UNIT = 12,FILE= ' CIRCLE / ,STATUS= 'UNKNOWN' ) 

OPEN (UNIT=13,FILE= ' CONTROL ' , STATUS= ' UNKNOWN ' ) 

READ( 1 1 , X ) M 
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READ(ll,X)f1 

READ(il,X)N 

READ(11,X)T 

READ(11,X)L 

READ<il,X)(UAL(I ),I=1,M) 

READ(11,X)< (UEC(I, J ), J = 1 , PI ), I = 1,N > 
READ( 13, X ) (XflASSC I),I=1,N) 

READ( 13, X )KTinE 
READ( 13,X )ALFA 
DO 1 1=1, N 
DO 1 J = 1,N 
C(I,J)=0 
D(I,J)=0 

IF < I . EQ ♦ J )C ( I , I ) = 2 ♦ XALFAXXMASS ( I ) 

IF ( I ♦ EQ . J )D ( I , I ) = ALFAXX2XXP1ASS ( I ) 

1 CONTINUE 
DO 2 1 = 1, PI 

0P1N= (0« ,-l ♦ )XUAL( I ) 

om(d=(0.,i. )xonn 

CEN( I ) = -ALFA+0f1( I ) 

LAH( I ) = CEN( I ) 

2 CONTINUE 

DO 7 IR = l,f1 
RAD( IR ) = 0 
DO 6 I S = 1 , H 
GRS = 0 

DO 5 J= 1 , N 
GJS = 0 

DO 4 1=1, N 

GJI = -(C( J,I )*LAPI(IS)+D(J,I ) )/Of1( IR ) 
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GJI=-(C(J,I )*LAM(IS)+D(J,I ) )/0M(IR) 

GJS=GJS+GJIXUEC( I, IS ) 

4 CONTINUE 
GRS=GRS+UEC ( J, IR )XGJS 

5 CONTINUE 

IF < IR . EQ . IS )GRS=GRS+ ( 2 . *LAM ( IR )XAIFA+ALFAX*2 )/0M ( IR ) 
IF(IR*EQ.IS )CEN ( IR )=CEN( IR )+GRSX0. 5 
IF ( IR .NE ♦ IS )RAD ( IR )=RAD( IR HCDABS ( GRS ) 

6 CONTINUE 

7 CONTINUE 
URITE(12,100)(LAM(I),I=1,M) 

100 F0RMAT(2X,' IDEAL EIGENUALUES'//,25(2X,2E15.5/ ) ) 

WRITE ( 12, 150 ) (XMASSC I ), I = 1,N ) 

150 F0RMAT(2X, ‘'REGIONAL MASSES V/,25(E15.5/) ) 

WRITE(12,200) 

URITE( 12,300 )(CEN( I ),RAD(I ), 1=1 ,M) 

200 F0RMAT(2X, 'NEIGHBOURHOODS OF THE CLOSED-LOOP' 

1 , ' EIGENUALUES ' /, /2X, 4X, ' CENTERS! FIRST-ORDER APPROX)', 5X 

1 2X 'RADII ' / ) 

300 F0RMAT(2X,2E15*5,5X,E15*5 ) 

CLOSE ( 11 ) 

CL0SEC12) 

STOP 

END 

BOTTOM 
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